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VAPORIZATION  0?  SOLIDS  BY  ABSORBED  RADIATION  EN3RQT 


E.M.  Shakhov 
(Moscow) 

Inzhenernyi  zhurnal  (Journal  of  Engineering)  1,  Mo.  1)., 

27-38  (1961) 

The  heating  up  of  a  solid  body  absorbing  radiant  energy 
was  discussed  in  earlier  communications  (Refs.  1,  2). 

When  the  density  of  the  radiation  flux  is  sufficiently 
high,  the  body  surface  may  acquire  a  very  high  tempera* 
ture  at  which  the  substance  will  inevitably  undergo  some 
physicochemical  transformations  (molting,  evaporation, 
decomposition,  etc.).  Let  us  consider  the  case  where  the 
body  absorbs  a  high-power  radiation  pulse  Q  of  short 
duration  7^.  By  the  <  nd  of  the  time  interval  7*  the  surface 
temperature  will  have  reached  a  value  proportional  to 
Q /Vti  when  Q  is  sufficiently  large  and  T  is  suf fi¬ 

de  ntly  small,  the  temperature  may  rise  to  any  conceivable 
level.  Therefore,  to  evaluate  the  heating  offect  of  powerful 
radiation  fluxes  abj  cubed  by  a  solid  body,  it  is  essential 
to  co  l dor  the  attendant  physicochemical  transformations  of 
the  substance.  Of  all  such  phenomena,  molting  is  the  simplest 
and  most  thoroughly  oxplorod  (Refs.  3,  li)  •  The  present  study 
is  conco  d  with  the  unidl  ensional  problem  which  Involves 
the  v  ^ori/.i  .ion  of  a  ;olid  body  h  \tid  by  the  lncldo.it  radiant 
«n  i.0y  It  f-i-.bs.  Ti  e  hbaor/.lon  Is  a  ;  iv  \e-l  to  O'-jor  on  the 


r  - ///i  »-2 


surface,  and  tha  radiation  flux  is  taken  to  bo  tho  equi¬ 
valent  of  a  heat  flow  q  along  the  boundary  whose  magnitude 
depends  on  surface  temperature  and  time  t.  Two  eases  of 

a 

evaporation  are  considered.  In  the  first  case  the  body. 
Initially,  is  contiguous  to  a  vacuum;  in  the  second,  the 
ambient  3pace,  at  some  distance  from  the  body.  Is  filled 
with  gas.  Both  vapor  and  gas  are  treated  as  perfect  gases. 
The  mechanism  of  vaporization  is  however  substantially 
different  in  the  two  cases  (Refs.  5,  6) .  When  evaporation 
is  intense,  in  the  case  of  a  body  initially  contiguous  to 
a  vacuum,  the  rate  at  which  the  vaporization  products  are 
removed  is  determined  by  the  velocity  of  expansion.  The 
latter  is  a  gas-dynamic  process,  which  is  sufficiently 
rapid  and  cannot  limit  the  speed  of  the  entire  procese. 

For  this  reason  the  loss  of  solid  mass  a  per  unit  of 
evaporating  surface  proceeds  at  about  the  same  rat#  ee 
evaporation  Into  a  vacuum  In  the  second  case,  however, 
where  the  ambient  spaoe  is  filled  with  a  gas-vapor  mixture, 
the  evaporation  products  are  removed  from  the  surface  by 
diffusion,  at  a  rats  which  is  low  as  compared  with  the  rat# 
of  evaporation  into  a  vacuum  -  the  highest  possible  for  a 
given  temper sture  of  the  body  surface.  7 hat  happens,  appa¬ 
rently,  is  that  most  molecules  rhich  have  left  the  body  end 
entered  the  gas  phase  return  to  "he  surface  and  once  again 
become  densely  packed  (condensation),  while  only  a  few  are 
permanently  reaoved.  It  follows  that  a  near  -  equilibrium  is 
Tfialntalr jd  in  the  pro-,  is,  and  the  state  of  vapor  envelop¬ 
ing  the  aurftce  is  close  to  saturation* 
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I •  Evaporation  fro»  the  surface  of  a  body  Initially  contlguout 

to  a  vacuum 

» 

1.  Vfe  shall  consider  the  relatively  simple  ease  o f  a  solid  — 
a  half-3pac«  x  <[  0  --  whose  Initial  temperature  Is  uniform  through¬ 
out.  The  body  absorbs  radiant  heat  along  the  boundary  with  an  inten¬ 
sity  q  dependent  on  surface  temperature  and  time.  Since  the  body 
occupies  a  negative  region  of  space,  the  heat  flux  q  should  also  be 
regarded  as  negative.  Let  0  denote  the  difference  between  the  Ini¬ 
tial  body  temperature  and  instantaneous  temperature  at  any  arbitrary 
point  x  of  the  body.  It  is  assumed  that  heat  propagation  within  this 
solid  obeys  the  linear  equation  of  thermal  conduction 


At  df*  ' 


(1.1) 


and  that  vapor  motion  is  described  by  equations  of  gas  dynamics.  On 
evaporation  the  body  surface  changes  lt3  position  in  specs.  Let  va 
designate  the  law  governing  its  motion  by  x0(t);  the  translation 
velocity,  by  x0(t);  and  by  J»0,  pQ,  TQ  and  uQ,  the  density,  pressure, 
absolute  temperature  and  ma3S  volocity  of  a  vapor  pa.  tlels  in  tho 
neighborhood  of  the  evaporating  surface.  ?roa»  the  laws  of  msss  and 
energy  conservation  the  following  expressions  con  now  bs  dorlved 
Interrelating  the  vn-cncrrn  parr-n*ors  of  the  evaporating  surface X 


—  py,  (0  -  {»» —  (0)  r=  *"• 


(1.2) 


(1.3> 

of 

Here  are  tho  density  end  coefficient  of  heat  conductivity 

of  a  solid  body,  and  t»  is  the  specific  heat  of  VAporixatlon  taken  t9 
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be  a  known  function  of  temperature. 

What  we  know  of  conditions  expressed  by  (1.2),  (1.3),  in  con¬ 
junction  with  initial  data  and  boundary  conditions  at  infinity,  is 
not  enough  to  determine  a  unique  solution  of  eq.  (1.1)  or  a  system 
of  gas-dynamic  equations,  because  our  knowledge  of  boundary  condit¬ 
ions  existing  on  the  evaporating  surface  x  m  xQ(t)  is  inadequate 
(Ref.  7).  Prom  the  standpoint  of  physics,  this  is  quite  obvious, 
since  conditions  (1.2)  and  (1.3)  not  reflect  the  specificity  of 
the  process.  Por  additional  relationships  governing  an  evaporating 
surface  we  must  turn  to  thermodynamics,  or  to  the  molecular-kinetie 
theory  of  the  structure  of  matter. 

We  shall  attempt  to  determine  in  what  way  the  mass  evaporation 
rate  depends  on  the  temperature  of  body  surface  and  the  state  of 
vapor  in  the  neighborhood  of  this  surface.  The  rate  of  any  physico- 
chemioal  process  can  be  expressed  with  sufficient  accuracy  as  the 
difference  between  the  respective  rates  of  the  direct  and  rewerso 
processes  (Ref.  8).  Por  the  case  of  evaporation,  we  haws 

m  -«  m,  —  m_,  (i.M 

where  4^  is  the  direct  process  rate,  or  the  rate  of  evaporation  into 
a  vacuum,  which  depends  on  surface  temperature  alone}  *.  is  the  rato 
of  a  reverse  process,  i.e.,  condensation,  dependent  on  the  state  of 
vapor  noar  the  body  surface.  The  values  of  and  can  be  calcula¬ 
ted  on  the  basis  of  statistical  physics.  Por  the  rate  cf  mass  evapo¬ 
ration  into  a  vacuum  we  have  (Ref.  9) 

m.  (I  ' r)  y&s • 


(1.5) 


where  T.^  is  the  absolute  temperature  of  body  surfaoe;  p#  i*  the 

saturated  vapor ■  pressure  corresponding  to  this  temperature;  R  is  the 
» 

g as  constant  cf  vapor;  r  Is  reflection  factor,  i.e.,  the  fraction  of 
the  total  number  of  molecules  that  reaches  the  body  surface  and 
undergoes  condensation.  To  calculate  ra_  we  3hall  assume,  as  did 
earlier  investigators  (Ref.  10),  that  an  equilibrium  distribution 
of  molecular  velocities  exists  near  the  body  surface,  the  average 
velocity  being  oqual  to  UqS 


here. 


—  w  • 

m.  ~  —  (1  —  r)  ^  ^  dv,  dv„  ^  fvM  dot, 


i  „ _ jl . exp  L  tv*~  l 

(2 t  2  RT,  ) 


Is  the  ..axwellian  function  of  velocity  distribution  .  Evaluation  of 


the  integral  yields 


m_  -  (1  —r) 

V'snKT, 


There  «,  -  -  Ai, i is  the  .'‘ach  number  near  the  body 

i  -XT,  Y  - 

surface;  f  is  the  ratio  of  the  spocific  boats  of  the  vapor; 

OS 

icrfc  x  -  -  —  -r*'-  x  ---  C  t’’  dx. 

v»  »'*  i 

For  the  mass  evaporation  rate  .a  now  obtain  the  expression 


in  m,  --  (l  -  -  r) 


yxAkft 


Y a  ierfc  v- 


(1.6) 


Frv-i  co-mltion  (1..1)  it  follows  that  since  -&■<!; 

-».(')  ■  * 

we  c  \n  th.  iroforo  write,  with  high  accuracy 


m  ■  ,vi#. 


(1.2‘) 
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Comparison  of  (1.2*)  with  (1.6)  gives 


m  = 


2u, 


2u»  +  ierfc  u. 


'=“»»+• 


(1.7) 


Consequently,  the  ratio  for  the  loss  of  solid  mass  in  the  two 
cases  -  evaporation  into  an  ambient  space  filled  with  pure  vapor 
as  against  evaporation  into  a  vacuum  -  will  depend  solely  on  the 
Mach  number  value  ;  oar  the  evaporating  surface.  Function  m/d4 
rises  monotonically,  and  very  rapidly,  with  the  rise  of  the  Mach 
number.  For  a  moncatomic  gas  ()f  =  5/3) »  when  Mq  »  0.5,  m/m*  m  0.805 
while  when  MQ  .  1,  it/m^  s  0 .96l|.. 

It  will  next  be  shown  that  the  Mach  number  cannot  be  much  less 
than  unity.  In  fact,  for  M0  ^  1  we  have  from  (1.4)  and  (1.7) »  with 
an  accuracy  limited  by  terms  containing  M0  in  the  zero  power  — > 

*+  —  a.;  i.e.,  the  state  of  vapor  near  body  surface  is  close  to 
saturation.  Let  the  parameters -of  saturated  vapor  be  denoted  by  the 
subscript  Evaluating  the  ratio  of  the  velocity  at  which  the  eva¬ 
poration  front  is  propagated  to  the  speed  of  sound  in  saturated 
vapor  will  give 


»>(!)  m 

0  rxa 

•  • 


_L_  A 

K&u  p 


This  means  that  whon  the  vapor  Is  in  a  nearly  saturated  state, 
the  front  of  an  expansion  wave  is  propagated  at  a  speed  which  sub¬ 
stantially  exceeds  the  velocity  of  a  shifting  evaporation  front. 
The  rarefaction  waves  will  reach  the  evaporating  surface,  lowering 
the  3poed  of  sound  and  simultaneously  reducing  the  gas  density 
near  the  surface.  This  in  turn  a ill  reduce  the  condensation  rate 
m_,  hence  both  mass  velocity  -.1  the  "ach  nu.bar  will  increase. 

It  follows  tv  it  should  the  .  ipor  re.iuh  a  3tate  close  to  saturation 
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(as  In  the  cate  where  the  vapor  region  is  Initially  limited),  this 
condition  will  soon  be  upset.  If,  on  the  other  hand,  the  vapor 
region  is  consistently  bound  by  a  vacuum,  the  vapor  will  nevar  reach 
a  state  close  to  saturation. 

Qualitative  cons idorations,  along  the  same  lines  .as  those  set 
forth  above,  suggest  that  in  general,  near  the  body  surface,  the 
velocity  at  which  the  front  of  an  expansion  wave  propagates  cannot 
materially  exceed  tho  velocity  of  a  3liifting  evaporation  front.  Sow, 
since 

-Uj  _  p -  (1  -  M,). 

-■>.(0  «*• 

it  follows  that  for  <T  1 

P  -i,(0 


Tha  Mach  nu  'bar  is  therefore.  siithor  greater  than  unity  or  (for 
Mo<l)  is  clora  to  unity.  Consequently,  on  the  basis  of  (1.7)  we 
can  approximate  4  .  This  assumption  scorns  all  the  more  justified 

since  the  reflection  factor  r,  which  enters  into  m^,  is.  determined 
experimentally,  with,  m  accuracy  not  exceeding  .*&%  (Ref.  11). 

In  .".imlng  up,  no  conclude  tbit  ‘he  tho r.nal  problem  on  tho  heat¬ 
ing  up  of  a  body  able  to  ovnyorata  can  bo  solved  independently  of 
determining  the  ms  motion,  and  is  formulated  as  follows.  It  is 
required  to  find  a  Uirnnud.  j  field  within  the  body,  by  way  of 
solving  tha  heat  cm  uctlvlty  equation  (1.1),  under  c  editions 
exist!  nr  along  an  un'cr.ovn  boundary  »0(t) 


<lf 


VlU 

L  ■ 
m 

m 


(1.8) 

(1.9) 
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and  under  the  initial  condition 

6(x.0)  =  0.  (1.10) 

Hero  is  the  value  of  function  9.g(x,t)  on  the  surface. -A  dis¬ 
tinctive  feature  of  this  nonlinearVproblem  stems  from  the  physical 
fact  that  a  substance  is  capable  of  evaporatirg  under  any  condit¬ 
ions  but  a  different  rate,  which  depends  largely  on  the  surface 
temperature.  The  peculiarity  here  is  that  this  surface  "temperature * 
0*(t)  is  not  preassigned  but  must  be  determined,  along  with  the 
entire  temperature  field  and  the  uknown  mobile  boundary. 

2.  We  shall  next  show  that  a  unique  solution  is  actually  possible 
for  the  problem  as  formulated  above.  To  this  end  a  solution  of  (1.1) 
(Ref.  12)  in  the  form  of  an  integral  is  given.  Under  conditions 
expressed  in  (1.8) -(1 .10) ,  it  is  written  In  the  fom 


M*).  *)♦(*;  •■)*»  + 

|  % 
+  \  [C<*.  t,  *,(*>.  v)  r.(t)  _  ***  (x.  I;  *,(*). 


(1.11) 


where 


Gfat,  fcx)- 


is  a  Emotion  expressing  the  effect  of  an  Instantaneous  point  source 
for  the  equation  of  heat  conductance  (1.1).  The  solution  of  (1.1)  ss 
given  by  (1.11)  will  be  unique  if  9#(t)  is  specified  for  the  boun¬ 
der/  as  s  slngle-vnl;\od  function  of  t  i<l  if  the  corresponding 
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integrals  actually  exist.  It  follows  that  8.~(t)  is  not  necessarily 
conjugated  with  initial  :i  Lions . 

Let  us  consider  too  intx-grals  which  appear  on  the  right-hand 
3ide  of  (1.11).  The  Integral 

I 

t;  xtKx),  t)-;l.(t;  0.)dr 

is  continuous  throughout,  including  the  line  xQ(t),  if  the  function 
q  is  differentiable  at  all  points,  except  perhaps  the  point  am  0, 
since  in  the  neighborhood  of  this  point  It  satisfies  the  inequality 

!  -7 !  <  const  t 

The  second  integral  In  the  right-hand  portion  of  (1 .11) 


in*,  o  (  \G(x.t,  «.(*).  0  f.(M  -  "•  -(r’  -LL*-— ---1  0. (t)  dx, 

;  L  di  J 

become  ill.  :mc»U3  as  ! itror.jo:  Vno  line  x  —  xQ  (t).  Here  the 

following  ■■  nation  will  '.v  •  d  (lef.  13)  for  any  continuous  function 
Qjt): 

t'(r.(0  0..')  V  (-,(0,0  ’-9.(0- 

"o  shall  nea  co.niL  .•  L  1'  it  for  x  x,(/)  -0  in  (1.11) .  At  the 
wo  i.  ve 

o. (0  x $  (V  v (0.  t;  rjt\  t,  ;  lf( r. (0.  0  +  “ ••  <0- 


Tf  a  s  ’■it*  -n  of  V-o  i  ■<  ists,  fan  function  0_(t)  relat¬ 


ing  t)  this  j  '.j;'-'.-,  *  -  1 1  ,  i  .!  f  •  j.mtion 


r 

(0  <■  •*(»).  ♦)••?(*;  ®J«lT  r 

;-  2^G(»,,0.0  ')-.(-)  *  je.(t)rft. 


(1.12) 


(1 .13) 
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Slnc«  G(xt(0,i;  Jt,(r).T)  and  ftllLll  at  *>*-.  t  aequira  a  singula- 

rity  of  tha  fora  —  ,  it  foil  ova  that  whan  q(t;8»)  haa  tha  font 

*  r<-i 


-p1-*.  r»g »  *hara  g(t;84)  ia  everywhere  a  continuous  function  of 
its  arguments,  tha  system  (1.12),  (1.13)  haa  a  uni qua  continuous 
solution.  For  a  sufficiently  small  value  of  t,  this  solution  can  ba 
‘found  by  tha  method  of  successive  approximations .  It  is  concluded 
directly  that  a  solution  of  (l.l)  exists  and  is  unique,  under  eon* 
ditions  specified  by  (1.8) -(1.10) . 

The  sane  thing  can  be  demonstrated  for  a  body  of  finite  thick¬ 
ness  .  The  value  0  in  this  case  must  be  replaced  by  the  Green's 
function  corresponding  to  the  boundary  conditions  on  the  underside 
of  a  plate. 

3.  We  shall  now  consider  an  exact  solution  of  the  problem,  for 
an  elementary  ease.  Let  us  solve  a  reverse  problem,  assuming  that  a 
constant  temperature  0*  .  const,  is  sstablished  instantaneously  on 
the  body  surface.  This  corresponds  to  a  heat  flow  q,  which  will  be 
determined.  When  the  temperature  of  a  surface  is  constant,  the 
velocity  of  its  translation  in  a  space  x^(t)  is  likewise  a  known 
constant.  Substituting  the  appropriate  values  of  8*  and  Xp(t)  into 
(1.12),  we  obtain 


&s.W+**0l 


t 

VT* 2HF5J 


(1.14) 


From  (1.14),  by  use  of  operational  calculus,  a  function  ^(t;0#) 
ia  derived,  characterising  a  heat  flow  pens trs ting  tha  body! 


"-at-  tl•lS, 


10 


The  magnitude  of  the  heat  flow  q  can  now  be  determined  from 
(1.9). 

In  the  ca'se  where  no  evaporation  takes  place,  xt (t)  r=  6 «.  0  ,  and  W£ 

obtain  the  familiar  formula  based  on  the  theory  of  heat  conductance 

9  a  -is-  >  which  In  our  case  holds  for  sufficiently  small  values 
YKS 

t.  When  t  — »  oo,  a  constant  heat  flow  directed  into  the  body  becom¬ 
es  stabilized: 


9(0o)“s 


*.(*>  , 

* 


(1.16) 


It  will  be  noted  that  the  same  formula  can  be  derived  directly 
from  the  equation  of  heat  conductance  (1.1),  if  the  solution  is 
sought  in  the  form  of.  a  heat  wave  propagating  at  a  uniform  speed 
«(jt,  0 «•  #(y)  ,  where  y  —  —  ~p-{x~-xJit)l). 

We  then  have 


Hence, 

J*  _ |  »«(Q  j 

a1*  T' 

For  y  OU^  i#(0*  this  expression  is  Identical  with  (1.16). 

Thus,  for  the  case  where  t  _ *  oo,  the  particular  solution  under 

consideration  goes  over  asymptotically  into  a  solution  based  on  a 
heat  wave  propagating  at  uniform  speed.  As  seen  from  (1.15),  the 
rate  at  which  this  thei-nal  process  approaches  s  steady  stats  is 
determined  by  the  ma.'ni iude  of  parameter 

r  2Vn  * 

The  velocity  of  the  shifting  evaporation  front  is  low,  aa  a 
rule,  and  the  thermal  diffuaivity  X  is  of  the  order  of  unity  (in 
the  CCS  system).  Therefore  the  velue  of  the  parameter  0T”  is  smell. 
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However,  since  the  tine  required  for  the  establishment  of  a  quasi* 
stationary  process  such  as  a  heat  wave  is  of  the  order  l/b,  the 
tine  period  involved  may  be  quite  large  (by  comparison  with  the 
characteristic  tlno). 

II.  Evaporation  into  an  ambient  space  filled  with  gas 
at  some  distance  from  the  body 

1.  Let  a  solid  half-space,  as  in  Section  I,  occupy  initially  a 
region  x  <0,  while  the  region  j>0  is  filled  with  a  gas.  It  is 
further  assumed  that  evaporation  of  the  substance  from  the  body 
surface,  due  to  absorption  of  radiant  energy,  sets  in,  beginnlx^ 
at  t  m  0.  The  state  of  a  gas  mixture,  it  will  be  remembered,  la 
determined  principally  by  the  processes  of  molecule  transfer  (sueh 
as  visoosity,  heat  conductance  and  diffusion).  In  accordance  with 
this  basic  concept,  we  shall  neglect  the  derivatives  of  pressure  in 
the  Havier-Stokes  and  energy  equations .  This  is  a  common  assumption 
in  the  theory  of  oc  lbustlom  (both  steady  state  end  non-steady  state 
processes)  (Ref.  14).  Unlike  the  processes  considered  in  this  work, 
combustion  involves  e  chemical  reset  ion.  Ve  shall  further  neglect 
both  thermal  and  baric  diffusion,  as  well  as  the  heat  flux  end  its 
dissipation  resulting  from  diffusion.  The  gas  peraaeters  can  new  be 
determined  from  e  system  of  differential  equations 


PJT +9mZT 


J. +£-* 

".(M  hi* 


(2.1) 


•12< 


where  e  Is  the  mass  concentration  of  vapor.  System  (2*1)  must  be 


solved  Jointly  with  the  equation  for  beat  propagation  in  a  solid 


body 


(2.2) 


under  the  specific  conditions  ormass,  concentration  and  energy 
conservation.  The  solution  is  derived  for  the  resulting  relation* 
ships  on  the  evaporating  surface: 

X  =•  x ,(/).  m  -  ?,;.(/)  -  P.K-x.(0)s=?^* 

(23) 

p,  (i  -  €.)  («,  -  i,  m  -  hd  ii^s!  _  & 

Subscript  0  refers  to  gas  parameters  near  the  surface;  subscript 
1,  to  values  characterizing  the  solid  body.  Conditions  on  the  eva¬ 
porating  surface  x  m  x0(t)  are  expressed  by  (2.3).  To  this  systea 
two  more  equations  must  be  added,  which  describe  a  thermodynamic 
equilibrium  between  the  solid  body  and  its  vapor: 

*-=*•(0.  Tj-r*  C.  -c.trjb  (2.1;) 

where  ea  is  the  coneentratlon  of  saturated  vapor  at  a  tempers  tors 
T0<  At  infinity  we  shall  prescribe  the  following  conditions: 

**-  +  «•.  « » 0,  c - 0,  e *-■  N*  T-f» 

r,-r^  (2.5) 


T.m  initial  conditions  will  bs: 

#~d ,  M -«(*,<& /*»«,  >-mjc,ov  r-rts,0i  r,-r,u.o>.  (2.6) 

In  place  of  x  s  new  Independent  variable  -  the  Lagrange  coordi¬ 
nate  -  is  lntrodueed:  n  *  -  «W  +  i'4* 
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1.  Simultaneous  determination  of  temperatures  and  coneentrationa 


n>-*(0 

0 

n<  — «(0 

* !-—«(/) 


ST 

« 


dtj*’ 


.  Sr 

.  ft  £Tj 
*  dtf  9 

TW* 


U<0  .-s-x.p.^'+Xp^-. 


1-e,  * 

n  -  ±.eo,  T  Tm.  e  =■  0, 

i-o.  r-r(jc,o).  e  =  o. 


r.-Tn. 

iV-r^o*. 


II.  Determination  of  velocity  and  density 


Sn-  SNl _ LA*  t 

.  si  **'  a*  si  “  * i* 

*«-*(/)  I 

n  »  a>.  e  - 0.  P"N» 

1-0.  u  «(x.O).  P  »P(*.  0*. 

The  subsequent  discussion  will  be  concerned  in  the  main  with 
the  fields  of  temperatures  and  concentrations. 

2.  To  cite  one  example  of  an  exaet  solution  of  the  problem*  let 
us  assume  that  at  the  initial  moment  a  half- space  x>0  Is  filled 
with  gas  whose  temperature  Is  uniform  throughout.  The  temperature 
of  the  solid  body  Is  likewise  uniform  everywhere! 

r(*.0)-r»  r,c*.o)-r,» 

Starting  from  t  a  0  the  solid  body  keeps  evaporating*  either 
because  It  is  as ted  upon  by  a  radiation  flux  equivalent  to  a  heat 
flow  of  the  font 


« 


(Q^O  possibly  depends  on  surface  temperature) ,  or  because  of 
mere  Interaction  with  the  gas,  since  the  body  and  the  gas  had  diffe¬ 
rent  Initial  temperature. 

From  considerations  of  similarity  it  is  concluded  that  the 
solution  of  problem  I  must  depend  on  a  single  dlmen3lonl5jpeS:$uarl» 
able  *  —  t 

r  = Tm+(rm-Tm)i(2).  «-«(*  v" 

The  movement  of  the  boundary  is  proportional  to  /t: 

—  m  [ t )  =*  n#  (/)  =*  ztY*i  (**  <  0). 

Assuming  *0,  9Q,  e0  to  be  known,  we  shall  find  for  functions 
0(z),  0^(z),  c(z)  expressions  which  will  satisfy  the  original  equat¬ 
ions,  the  initial  conditions  and  the  conditions  of  infinity! 


trie 


•'«- !  -P-M* 


2 -eric 


1£55. 


trk 


2VJH 


c(z)- 


2 -eric 


aVTQw 


(2.7) 


eric 


a/isT 


•*k 


where 


crfc*--i-C  s-**d*. 

/»  i 

Since  for  »  it  follows  that  *-*(*).  • 

remaining  unknown  values  z0,  0Q  are  determined  from  relationships 
expressing  the  laws  of  energy  and  component  conservation.  These,  In 


ou~  case,  have  the  font 

-  kc  -  -9<w-*;w +*•  Wb 


(2.8) 

(2*9) 
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wh  ere 


2X,p,  (r.-r,.) 


.  S« 


L#  is  ths  hsat  of  evaporation  at  Ta  *joo* 

We  shall  now  transform  the  system  (2.8),.  (2.9),  introducing 
notations 


”w— '£l"(*,fcI^)-VS 


•(-&) 


-*-3= 

*YfM 


>flL 


»»W 

We  than  haws 


Considering  that  the  specific  heat  of  the  body  c^  and  that  of  the 
vapor  0^  are  constant,  the  dependenoe  of  the  heat  of  evaporation  on 
temperature  will  be  linear! 

— *a)(7»— Tia). 

Consequently, 


t<»4  -*  l  + 


(2.10) 

By  use  of  (2.10),  equation  (2.8)  can  be  rewritten  in  the  font 


KM 


adtjaeadt 


•<*) + w- *  V-S  (r.-r^>% 


(2.U) 


If  q  is  assumed  independent  of  surface  temperature  (0  *  const.), 
the  system  of  two  transcendental  ecuatians  (2.9),  (2.11)  will  be 


reduced  to  a  single  equation  which,  on  letting  #— 


,  can  be 
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be  written  aa 


y(l  +  etly ) 


J _ *• 


/r  i-«.  ‘ 

« 

since  eric (—y)  =•  1  +  erf  y.  .  This  equation  has  a  solution,  which  la 
unique. 

*7e  can  now  solve  problem  71,  which  involves  the  distribution  of 
velocity  u(r),  t)  and  of  density.  A  solution  is  found  In  the  font 


Calculation  yields 


•  *»~** 

"  (1.  0  =*  *  *  » 

p  • 


where  c(z0)ls  determined  by  substituting  uQ  and  fQ  Into  the  condit¬ 
ion  of  mass  conservation  (plut^m)  .  We  obtain 


The  problem  under  consideration,  it  will  be  noted,  remains  self- 
similar  also  for  the  case  of  an  arbitrary  dependence  of  ®  on 

tempera ture  and  concentration. 

3.  Problem  I  is  a  relatively  complex  nonlinear  problem,  hence  an 
exact  solution  of  it  can  apparently  be  found  only  for  some  particu¬ 
lar  cases.  The  above  self-similar  problem  la  the  simplest  of  these 
eases.  When  the  equations  are  of  the  same  type  as  the  equation  for 
heat  conductance.  It  is  often  possible  to  obtain  another,  quit# 
simple  yet  rather  interesting,  solution,  based  on  a  wave  propagating 
at  uniform  speed.  We  cannot  however  consider  here  the  ease  of  a  wave 

-18- 


propagating  In  one  direction  only,  since  a  solution  of  equations  In 
this  form  will  not  satisfy  all  the  boundary  conditions  involved. 
Inasmuch  as  deriving  an  exact  solution  under  arbitrary  initial  and 
boundary  conditions  presents  a  highly  complicated  problem,  we  pro¬ 
pose  ^a^m.ethod  for  arriving  at  an  approximate  solution.  We  shall 
transpose  the  boundary  conditions  for  the  case  of  rj  m  to  the 

line  —  0.  Froblem  I,  considerably  simplified  thanks  to  this  assum¬ 
ption,  can  be  formulated  as  shown  below: 


0.<n<oo. 

—  oo  <  n  <o. 
'1=0,  r,=r* 


9T 

m 


it  -9* 

**  #*’  • 

uBk  ' 
h*e'  ■ 


—  L  • 


t  =  0. 


r,  =  r,(*.  0). 


The  oonditionm (0** - — -^-servos  to  determine  the  law  governing 

I—  d#  0% 


the  shifting  of  the  boundary,  after  the  fields  of  temperature  end 
concentrations  have  been  defined.  W«  are  now  ready  to  solve  our 
problem.  Let  ue  denote  by  (n. t; g, t)  the  Green  function  of  the  first 
boundary  prdblem  in  the  heat  conductance  equation,  where  /  Is  the 
coeffloient  et  the  second  derivative  with  respeet  to  for  tho  half- 
lino  *>  0s 


Functions  *•(«».<;  **)  are  obtained  from  &(%<;»*)  by 

substituting  fi  i  and  6  for  fi  .  Assuming  the  surface  temperature 
r(0,Q-r#W.  to  be  known,  we  shall  obtain  for  the  functions 
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fOl. 0,  e (»!,/),  Ti (n,0  expressions 


T  (1. 0  =  S  ^ <*•  0)  <n.  <;  y.  0) <*y  +  p  $  r, (x) 

c(r\.i)=o{ct 

S  d» 


*t%  (n.  /;  0,  t) 


(2.12) 


(2.13) 


r»(*t.<)~  ^r,(y.O)g>  (t|.  /;  y.  0) dy  —  ft,  ^  T, (t)  — <*X.  (2.1U)  J*t'’ 

■"her a  c,  =»c,(r.)  .  The  surface  temperature  T  (t)  is  determined  from 
the  condition  expressing  the  heat  balance  at  the  boundary,  fle  shall 
write  it  in  the  form  integrated  with  respect  to  time  t: 


(2.15) 


The  Integral  on  the  loft-hnnd  3ide  of  (2.15)  can  be  rewritten  in 


the  fom 


(2.16) 


where,  at  the  first  approximation. 


We  shall  now  set  up  the  partial  derivatives  of  functions 
1*(H.  A,  cfa 0,7*1  (n.t).  with  rospect  to  q.  .11020  functions  are  deter- 
tinod  from  eq.  (2.12) -(2 .1!;)  .  #e  shall  integrate  them  with  respect 
to  t,  and  shall  obtain  at  the  limit,  for  q  _*  Os 
"(n. x>  a\  _ f  c  <V°I  *  *  °>  '  r  , 


*  si  **  i  * 

,im  <£X“*  -U-r^— 

.  +$7,(0— 


.  St 


«**.  (it  t;  a  9 


d#A+  • 


(2.17) 

.  (2.18) 


(2.19) 
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Substituting  (2.17)-(2.19)  into  (2.15),  »•  obtain  by  use  of 
(2.16)  an  Abel  equation  expressing  the  function 


<HT.)  =  r.(0  +  Yr 

Wh  Yl'  . 


and  thus  arrive  at  the  following  equation  for  de 
temperature 


surface 


where 


f  4 1 


V*  <1  - 1) 


(2.20) 


„  •  »e  ,(0.>:jr.o) 

^  ~  r.)  -r  xlP,  J  r,  o,.  o>  »•  ^  dy-  • 

“  dg  (0.  <;  If.  0) 

-Xp^r(y.O) . - dy. 


If  the  value  of  the  heat  flux  does  not  depend  on  surface  tempe¬ 
rature,  then  F  is  a  function  of  tine  alone,  and  formula  (2.20) 
gives  an  implicit  expression  for  surface  temperature. 
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